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The conformal thin sandwich (CTS) equations are a set of four of the Einstein equations, which 
generalize the Laplace-Poisson equation of Newton's theory. We examine numerically solutions of 
the CTS equations describing perturbed Minkowski space, and find only one solution. However, 
we find two distinct solutions, one even containing a black hole, when the lapse is determined by 
a fifth elliptic equation through specification of the mean curvature. While the relationship of the 
two systems and their solutions is a fundamental property of general relativity, this fairly simple 
example of an elliptic system with non-unique solutions is also of broader interest. 
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PACS numbers: 04.20.Ex, 04.20.Cv, 04.25. Dm 

Within a space-plus-time decomposition 0, 0], Ein- 
stein's equations, just as Maxwell's equations, split into 
evolution equations and initial value equations. The lat- 
ter constrain the initial data and are generally solved 
by the use of potentials that satisfy elliptic equations. 
These potentials are expected to be unique, given suit- 
able boundary conditions. For the Einstein constraints, 
indeed, in all cases in which existence of solutions has 
been proved, the solution was also shown to be unique. 

However, we demonstrate in this paper that an impor- 
tant method of solving the Einstein constraints through 
potentials, the extended conformal thin sandwich (XCTS) 
equations ||, 0], admits non- unique solutions. Two so- 
lutions appear to exist even for arbitrarily small per- 
turbations away from Minkowski space; one solution is 
just perturbed Minkowski space, but the second one even 
contains a black hole for sufficiently small perturbation. 
The XCTS system thus provides a fairly simple example 
of a nonlinear elliptic system with non-unique solutions, 
which may be of interest outside general relativity. Be- 
cause the XCTS system lies at the heart of most modern 
schemes to construct binary neutron star [fj, |fj|, [21 la] and 
binary black hole |Tj, LLfJ, [ill LL2J initial data, our result is 
also relevant to the astrophysical problem of computing 
inspiral waveforms for compact objects. 

In general relativity, the initial data are the induced 
Riemannian metric gij and the extrinsic curvature K%j 
on a spacelike hypersurface. The data {gij, Kij) must 
satisfy four nonlinear constraint equations, and there are 
now two equivalent, complete approaches with which we 
can solve the constraints. One methods seeks (gij, Kij) 
directly while the other, the "conformal thin sand- 
wich" (CTS) method 0, seeks to construct (gij,d t gij). 
In general relativity, this difference is not trivial: Ki j de - 
pends only on the embedding of the slice (see, e.g. |F3j) 
while dt g^ depends also on the ambient coordinate neigh- 
borhood of the slice. Both methods result in four coupled 
elliptic equations, and many existence and uniqueness re- 
sults are known |14j . 

The extended conformal thin sandwich formalism adds 



a fifth elliptic equation for the lapse-function, which 
arises from specification of the time derivative of the 
mean curvature K = K^g 13 . In the extended system, 
the free data consists entirely of "variable & velocity" 
pairs (<?, q) 0], as desirable for a Lagrangian viewpoint, or 
a thin-sandwich viewpoint. However, the fifth equation 
couples strongly all five equations of the extended sys- 
tem, complicating mathematical analysis, so that to our 
knowledge, there are no rigorous mathematical results for 
the extended system, neither for existence nor for unique- 
ness. Similar difficulties, as well as the non-uniqueness 
results reported in this paper are not present for the stan- 
dard initial value problem, but arise when going beyond 
the mere initial value problem by placing demands on the 
slicing in the ambient spacetime via specification of dtK. 
Such issues connected with general relativity have been 
called traditionally the "problem of time." 
The CTS equations can be written as 

VV - - + ^-'A^A^ = 0, (1) 
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Here, V and R are the covariant derivative and the trace 
of the Ricci tensor associated with the conformal metric 
gij , which is related to the physical metric by g^ = ip 4 gij . 
Furthermore, 
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represents the conformal trace-free extrinsic curvature, 
(Lj3) i3 = 2V { - i (3 3) - 2/3 g i3 %(i k is the conformal longitu- 
dinal operator, and the traceless tensor tty =dtg%j repre- 
sents the time-derivative of the conformal metric. Once 
the free data (gij,Uij;K,N) are chosen, Eqs. and 
are elliptic and determine the conformal factor ijj and the 
shift-vector (3 l . The conformal lapse N is related to the 
physical lapse by N = -0 6 7V. 

The elliptic equation for the lapse, which follows from 
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the Einstein evolution equation for K, can be written as 

V 2 (AV) - (N^){^R+^K^+ 7 -^- s A^A l ,) = 

-4> 5 (d t K - f3 k d k K). (4) 

Viewing dtK instead of N as freely specifiable, the free 
data becomes {gij,Uij]K,dtK), and one must solve the 
coupled elliptic system Eqs. J2J) and (gj. This is the 
XCTS system. 

We will show non-uniqueness by explicitly construct- 
ing two solutions for a certain choice of free data. This 
choice is based on a linearized quadrupolar gravitational 
wave and follows very closely [lif. We set 



9ij — fij + A hij\ t=Q 



= A TF 5 d t h 



a u tiHj | t _ > 
K = 0, 

N = 1, for CTS 
d t K = 0, for XCTS 



(5a) 
(5b) 
(5c) 
(5d) 
(5e) 



where fij represents the Euclidean metric, TF§ means 
the trace- free part with respect to gij , and 

h i jdx i dx j = -8rK cos 8 sin 2 6drd0 - 2r 2 L sin 3 8d^d6. 

The functions K and L are are given by 

K = r~ 2 F {2) - 3r~ 3 F {1) + 3r~ 4 F, 
L = -r- 1 F {3) + 2r~ 2 F (2) - 3r~ 3 F (1) + 3r" 4 F, 



where F = F(r ■ 
profile, and F^ 



t) 



-(r+t-r ) 2 /w 2 denotes the radial 



d n F(x)/dx n . For small A, this choice 



of hij corresponds to a localized incoming gravitational 
wave at radius = 20 with width w = 1, having odd 
parity and azimuthal quantum number M = 0. 

Being based on a solution to the linearized Einstein 
equations, the free data JSJ satisfy the conformal thin 
sandwich equations QJ, and to linear order in 
A, but violate them to 0(A 2 ). Consequently, the XCTS 
equations have a solution which corrects the free data by 
0(A 2 ), as demonstrated in 0|. However, this is not the 
full story. 

We employ the spectral elliptic solver described in |J3 
to solve Eqs. JTJ, (0 and Q in a computational domain 
with outer radius 10 8 with Dirichlet boundary conditions, 
if> = N = = 0. For smooth problems, like the ones 
considered here, such a spectral method results in ex- 
ponential convergence, and allows construction of highly 
accurate solutions. All solutions presented in this paper 
converge indeed exponentially; generally, the maximum 
residual of both Hamiltonian and momentum constraint 
is < 10~ 6 , and other quantities we mention below have 
converged to a similar level. The convergence rate dimin- 
ishes along the "upper branch" of the extended system 
at amplitudes A < 0.2, but the relative errors remain 




FIG. 1: Solutions of the standard (four) conformal thin sand- 
wich equations vs. the conformal amplitude. As the critical 
amplitude A c ,4 ~ 2.01 is approached, the conformal factor 
grows to infinity. 



<10~ 4 even at A = 0.02. In all figures, numerical errors 
are much smaller than the line-thickness, so that essen- 
tially the "true" analytical solutions are plotted. 

Figure n presents solutions to the standard conformal 
thin sandwich equations for various amplitudes A. The 
conformal factor ip diverges as A approaches some crit- 
ical amplitude. For our particular choice N = 1, K = 0, 
Eqs. m and reduce to a set of equations which is 
very well understood mathematically |18l HIH: Solu- 
tions are unique and exist iff the Yamabe constant of 
is positive. Figure ^ suggests that the Yamabe constant 
of the metric Eq. I|5a(l is positive for small A and changes 
sign at A — Ac,4, similar to an example given in 20]. 

The situation is very different when solving the ex- 
tended system. Figure |21 presents the plot analogous 
to Fig. n for this case. As in the standard system, we 
have not been able to find solutions beyond some critical 
amplitude. However, here the similarity ends. For the 
extended system, the critical amplitude A c ,5 ~ 0.30422 
is much smaller than for the standard system. There- 
fore, existence of solutions cannot depend simply on the 
sign of the Yamabe constant of g^. Furthermore, as 
8 A = A c ,5 — A becomes small, max?/> varies as (5 A) 1 / 2 
and approaches a, finite limit ij) cr it ~ 1.0999 (the numeri- 
cal values for A c ,5 and ip C rit were computed by fitting a 
parabola). The log-log plot in the inset of Fig. ^confirms 
these statements. Inspection of the solutions shows that 
as (5^4 — all variables vary in proportion to (S A) 1 / 2 at 
any spatial coordinate location x, 

U-(A;x) « u crit (x) - v crit (x) \/~8A. (6) 
Here, u_, u cr n and v cr a represent the vector of all 
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FIG. 2: Solutions of the extended (five) conformal thin sand- 
wich equations vs. conformal amplitude. The conformal 
factor remains finite as the critical amplitude A c ,5 ~ 0.304 
is approached, as confirmed by the inset (SA = A c ,5 — A, 
Sip = tpcrit — max(VO). 



five variables (?/>, /3 l , iV), and the notation indicates that 
w_ depends on both A and the spatial coordinates x, 
whereas u cr u and v cr u are independent of A sufficiently 
close to A c ,5- 

Because of the parabolic behavior close to A C} 5, we 
expect a second branch of solutions, generalizing Eq. © 
to 



u±(A; x) w u crlt (x) ± v crlt (x) 



(7) 



The solutions found so far in Fig. [3] were obtained with 
trivial initial guess tp = N = 1, f3 l = 0, and consti- 
tute the lower branch it_ only; to converge to w + , an 
initial guess sufficiently close to u+ is necessary. Dif- 
ferentiation of Eq. (JSJ with respect to A yields v cr a = 
2{5A) 1 / 2 du_/d„4, so that, sufficiently close to -4 C ,5, 



-{Ax) 



XAx) 



ASA du - {A]x) 



dA 



(8) 



Approximating du-/dA by a finite-difference quotient of 
solutions with *4 = 0.300 and A — 0.301, we compute the 
right-hand-side of Eq. JSJ) and use it as initial guess for 
the elliptic solver at A — 0.300. The solver now converges 
to a different solution; for example, max?/; ss 1.12. This 
solution u+ for one A is then used as initial guess for 
solutions with neighboring A and repeating this process, 
we move along the upper branch and compute solutions 
u + for a wide range of A. Eventually, Fig. |3] emerges, 
which shows clearly the two distinct branches merging in 
a parabola at A c ,5- 

The solutions u + are quite remarkable. They exist 
over a wide range of A, down to small A. As A de- 



FIG. 3: Maximum of the conformal factor, minimum of the 
physical lapse N = ip 6 N , and maximal magnitude of the shift 
vs. conformal (unphysical) amplitude. Solid/dashed curves 
represent the lower/upper branch. The insets show enlarge- 
ments around the critical point; circles/diamonds denote in- 
dividual solutions along the lower/upper branch. 



creases, ip increases everywhere; the increase is particu- 
larly pronounced close to the maximum of tp which oc- 
curs in a ring in the equatorial plane at radius close to tq. 
Consequently, with decreasing A, the maxima of ip be- 
come increasingly pronounced and concentrated around 
r ~ ro in the equatorial plane. This necessitates higher 
angular resolution in latitude and diminishes the conver- 
gence rate of the spectral expansion, so that we choose 
to stop at ^4 = 0.02. We note, however, that we do not 
see any fundamental indication that the branch u + ter- 
minates; it appears that m+ extents to arbitrarily small 
A, with ever increasing -0. The physical (ADM) energy 
E of the computed initial data sets varies approximately 
proportional to max-0 — 1 along u+ (see also Fig. 0J. 
At A = 0.02, E = 50.012 (in units in which r = 20). 
This data set also contains an apparent horizon with 
mass M= y/ Area/ 167T = 49.872. The apparent horizon is 
oblate with equatorial radius 28.8 and polar radius 21.1. 
Apparent horizons are present for A ^ 0.025. 

Figure [3] shows clearly that the extended system ad- 
mits two solutions for the same free (conformal) data. 
However, the physical properties of the initial data sets 
depend on the physical quantities {g%j, Kij); for example, 
the metric takes the form 



(9) 



This can be interpreted as the superposition of a gravita- 
tional wave hij with physical amplitude A tp 4 on a physi- 
cal background tp 4 fij. The conformal factor ip is a func- 
tion of space, and therefore A "0 4 is also spatially depen- 
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FIG. 4: Energy E vs. unphysical amplitude A and physical 
amplitude .Amax. The left inset shows an enlargement around 
the critical point. The right inset plots the data of Fig. 0vs. 
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dent. To simplify, we take the maximum of ip, which 
coincides with the location of the gravitational wave at 
radius tq = 20, and define a representative physical am- 
plitude by .Amax = A max yb A . 

Figure ^ plots energy E both vs. A and Amax; the 
qualitative difference between both definitions of ampli- 
tude is apparent: The unphysical amplitude changes "di- 
rection" along the sequence of solutions leading to the 
largest energy at smallest A. In contrast, A max is con- 
tinuously increasing, precisely one solution corresponds 
to each value of A. m ax, and large energies are obtained 
at large A. m ax- The family of solutions passes smoothly 
through the critical point, and even for arbitrarily large 
Amax, a solution seems to exist. Thus, Amax represents 
the amplitude of the physical perturbation much more 
faithfully than A. 

We see that uniqueness of the problem is not a straight- 
forward issue, ft depends on the question: for a given set 
of free data (this is, for a given value of A), the ellip- 
tic equations may have two solutions (A < A c ,s) or no 
solution at all (A > A c ^). However, for each value of 
the physical amplitude A m ax as described, precisely one 
solution exists within the family of solutions considered 
here. 

While the relationship between uniqueness and non- 
uniqueness of the XCTS system as discussed above is in- 
triguing, the major result of this work lies in the fact that 
the XCTS system may have two solutions for the same 
free data. Besides the obvious importance of this result 
to workers in mathematical relativity, we note that this 



issue may complicate construction of astrophysical com- 
pact object initial data, which relies on some form of the 
XCTS equations H0QIEHHE1E3- Furthermore, 
the most common elliptic gauge conditions |2lj include 
precisely Eq. (QJ to determine the lapse-function. 
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